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In this paper we suggest a new approach for the multichannel Coulomb scattering problem. The 
Schrodinger equation for the problem is reformulated in the form of a set of inhomogeneous equations 
with a finite-range driving term. The boundary conditions at infinity for this set of equations have 
been proven to be purely outgoing waves. The formulation presented here is based on splitting 
the interaction potential into a finite range core part and a long range tail part. The conventional 
matching procedure coupled with the integral Lippmann-Schwinger equations technique are used 
in the formal theoretical basis of this approach. The reformulated scattering problem is suitable 
for application in the exterior complex scaling technique: the practical advantage is that after the 
complex scaling the problem is reduced to a boundary problem with zero boundary conditions. The 
Coulomb wave functions are used only at a single point: if this point is chosen to be at a sufficiently 
large distance, on using the asymptotic expansion of Coulomb functions, one may completely avoid 
the Coulomb functions in the calculations. The theoretical results are illustrated with numerical 
calculations for two models. 
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I. INTRODUCTION 

The three-body Coulomb scattering problem remains 
both a formally as well as computationally challenging 
problem An understanding of three-body scatter- 
ing has implications to many fields of science, for exam- 
ples, to combustion studies 0, reactions in the interstel- 
lar medium plasma chemistry [ij and other types of 
gas phase reactions [5[. Besides their fundamental in- 
terest, these applications have also generated significant 
experimental efforts. The double electrostatic storage 
ring (DESIREE) currently being constructed at Stock- 
holm university @ is designed to investigate the gas- 
phase reactions of oppositely charged species, of which 
three-body scattering reactions are the experimentally 
simplest subset, so generating high-quality data which 
will challenge theoretical studies as well as providing 
much-needed input to models of interstellar and plasma 
reactions. 

In Ref. we initiated a set of studies with the aim 
of obtaining a method for accurately computing state 
selective three-body multi channel scattering which also 
included Coulomb interaction. Our method was inspired 
by the mathematically sound approach of Nuttall and 
Cohen [8] related to exponentially decreasing or finite 
range potentials. Rescigno et al. Q modified this ap- 
proach to scattering by non-Coulombic but long-range 
potentials. The recent formalism developed in [7[ pro- 
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vides the mathematically solid basis for application of 
the complex scaling method [1, [t| to the single channel 
Coulomb scattering problem. 

In Ref. fiol ] we outline how this new formalism can be 
extended to problems in three-body scattering which in- 
clude Coulomb interactions. This extension is based on 
the same formal and numerical technique as was used for 
computing three-body resonances fill ] . Furthermore, this 
extension can be also combined with our three-body res- 
onance methods [ll|, E3 into a technique by which one 
may quantitatively identify the influence of resonances 
[13| on the three-body scattering cross section. The pur- 
pose of the present contribution is to study in detail the 
advantages and limitations of the present formulation of 
two-body multi channel Coulomb scattering. We are cur- 
rently investigating the possibility to extend this formu- 
lation to one that can be generalized to accurate studies 
of three-body multi channel scattering where collisions 
between charged fragments occur jg|. 

Formally, the new formalism is based on splitting the 
entire potential in a sum of two sharply cut-off poten- 
tials: a core potential, which is then of finite range; 
and a tail potential, where all but the Coulomb inter- 
action can be neglected. In a recent contribution [l4| 
we have studied the structure of the solutions of the 
three-dimensional Schrodinger equation for sharply cut- 
off Coulomb potentials. The derived formulation of the 
three-dimensional driven Schrodinger equation for prob- 
lems involving Coulomb interaction is shown to open the 
way for the forthcoming applications in three-body sys- 
tems as outlined in Ref. [lOj . 

The report is structured as follows. The main ideas be- 
hind this work are described in section Hi Al We begin by 
presenting the two-body multichannel scattering problem 
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in subsection lll Bl The analytically solvable equations for 
the diagonal part of the potential tail are investigated in 
subsection III CI while the formulation of the problem for 
the total potential tail is found in subsection III Dl The 
complex scaling theory is combined with the scattering 
theory in subsection lll El The numerical implementation 
of our theory is presented in section IIIII Two different 
examples are considered in subsections IIII Al and IIII Bl 
Finally, section HVl presents a summary of the report. 



II. THEORY 
A. Theoretical background 

The boundary conditions for a scattering problem are 
conventionally specified as the superposition of an in- 
cident wave and an outgoing wave. For two-body sys- 
tems with central interactions the Schrodinger equation 
can be reduced to a set of one-dimensional partial-wave 
equations. The boundary conditions for a single one- 
dimensional equation can easily be specified and numer- 
ically implemented. For the three body scattering prob- 
lem the boundary conditions at infinite distances between 
the particles are much more complicated. As such it is 
difficult to implement these boundary conditions in prac- 
tical calculations. 

One of the techniques commonly used to avoid the 
problems arising from these boundary conditions is the 
complex scaling transformation fl5| method. In this 
method one maps the radial coordinate r onto a path 
g a (r) in the upper half of the complex coordinate plane 

r ' — > 9a{r) , < a < ir/2 , 

9u{r) ~ const + r exp(ia) , as r —¥ oo . (1) 

This method is widely and successfully used when calcu- 
lating resonance energies in problems where the bound- 
ary condition at infinity is a purely outgoing wave. It is 
well known that the transformation (TT]) converts a purely 
outgoing wave to an exponentially decreasing function. 
Thus, if the scattering problem is reformulated in such a 
way that the boundary condition at infinity is a purely 
outgoing wave then, after a complex scaling transforma- 
tion, one obtains a problem with zero boundary condi- 
tions at infinity. 

This approach to scattering problems was first re- 
ported by Nuttall and Cohen in 1969 @. These authors 
suggest applying the Hamiltonian operator to the differ- 
ence between the scattering solution of the Schrodinger 
equation and an incident wave to yield an inhomogeneous 
(driven) Schrodinger equation. This difference between 
the functions apparently behaves as a purely outgoing 
wave at infinity. Hence, this equation is suitable for com- 
plex scaling, and @ used the uniform complex scaling 
approach 

9a(r) = rexp(ia) . (2) 



However, it appears that this approach is limited to the 
cases of finite-range and exponentially decreasing poten- 
tials. This is due to the following reason. The inhomoge- 
neous (driving) part of the driven Schrodinger equation 
is the product of the potential energy and the incident 
wave. The incident wave after the complex scaling trans- 
formation becomes a superposition of the increasing and 
decreasing exponential functions at infinity. Hence, the 
entire inhomogeneous part diverges. The class of finite- 
range and exponentially decreasing potentials form an 
exception where the method of Ref. [8] can be success- 
fully applied. 

The formulation of Nuttall and Cohen was modified by 
Rescigno et al Q for one-dimensional single-channel long- 
range potentials (explicitly except Coulomb potentials). 
These authors also start from the equation for the differ- 
ence between the scattering solution of the Schrodinger 
equation and an incident wave. However, instead of the 
potential V(r) they use the finite-range potential Vn(r) 
defined as 

Vr ^ = ( ,r>R. (3) 

The driven Schrodinger equation with this potential does 
not experience difficulties with divergence after apply- 
ing complex scaling since this potential is of finite-range. 
Furthermore, the solution of the unsealed problem with 
the truncated potential Vr(t) approaches the solution of 
the original problem with the entire potential V(r) as 
R —> oo. As shown in [j| the same is not true for the 
scaled equation. The solution of the scaled equation with 
the potential Vn(r) gives the incorrect scattering ampli- 
tude since the function is not analytic [lj| . Therefore, 
the authors of [9j suggested to use exterior complex scal- 
ing \Tt\ . This transformation belongs to the more 
general type given by Eq. (fT]). It is defined by the func- 
tion g a ,Q where Q is the scaling point. The inner interval 
[0, Q] is mapped onto itself 

9a,Q(r)=r, r<Q. (4) 

If Q > R, the solution of the scaled equation with the 
truncated potential Vfj(r) approaches the solution of the 
scaled equation with the entire potential V(r) as R — >• oo. 
One can calculate the scattering amplitude with desired 
accuracy by choosing a proper value of R. 

This method is not applicable to the Coulomb scatter- 
ing problem directly since as it is well known that trunca- 
tion of the Coulomb potential leads to noticeable errors 
for any truncation radius R. In the two body scattering 
problem the Coulomb potential can be implemented into 
the discussed approach if it is included in the free-motion 
Hamiltonian, while V(r) describes the short-range part 
of the interaction. In this case the incident wave is rep- 
resented by a Coulomb wave function, which is known 
analytically. This approach has been successfully used 
for calculations in atomic [l8| and nuclear [l9| physics. 
Unlike the two body case, an analytic solution for the 
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Coulomb problem does not exist if three or more parti- 
cles are involved in the scattering process. 

In our recent report Q we have shown how the method 
of exterior complex scaling can be generalized to the 
Coulomb scattering problem. Instead of truncating the 
potential we represent the entire potential as V(r) — 

(r) + V R (r), where Vr(t") is the same as in Eq. ^ 
and the potential tail V R (r) is given by 



V R (r) = 



, r<R , 
V(r) , r > R . 



(5) 



The approach discussed in this recent report 0] is based 
on solving the problem for the potential tail V R (r) at 
the first step. The solution of the scattering problem for 
the potential tail V R (r) plays the role of the incident 
wave. On subtracting this incident wave from the scat- 
tering wave function we obtain a function which asymp- 
totically behaves as a purely outgoing wave. After the 
transformation (JT|) with the function g a ,Q satisfying 
the boundary problem for this function has the trivial 
zero boundary conditions both at the origin and at infin- 
ity. 

In the present study we proceed with a formal as well as 
numerical study of the two-body single and multichannel 
problems. All the potentials discussed in this treatment 
that are denoted with a subscript or superscript R are 
defined analogously to Vr and V R in Eqs. §5§ and ([5]), 
respectively. 

B. The two-body multichannel scattering problem 

In the following discussion, consider the two-particle 
multichannel scattering problem with M channels. We 
assume that the interaction between the particles V nm 
(n, m = 1, . . . , M) depends only on the inter-particle dis- 
tance r and when r — > oo can be asymptotically repre- 
sented as 



V„ m (oo) 



(6) 



The quantities t n = V„ n (oo) are called the thresholds. 
The total interaction is given by the sum 



+ Km(r) + S n 



(7) 



The first diagonal term corresponds to the Coulomb in- 
teraction while V* m (r) describes the short-range inter- 
action which is assumed to decrease faster than r -2 for 
large particle separations. More precisely V nm (r) should 
obey the condition 



dr(l + r)|K ro (r)|<oo. 



(8) 



The partial wave multichannel Schrodinger e quat ion for a 
given angular momentum £ (see for example [2CH22j) has 
the form of a set of equations for partial wave functions 
*},(r) 



dr 2 



1(1+1) 2k f r, f 

T 2 T 
M 



- k 2 
fc / 



(9) 



71=1 



Here the coupling terms are of the form V nm (r) — 
(2fj,/h 2 )V nm (r), while the momentum k n and the Som- 
merfeld parameter r\ n in the n-th channel are defined 
through the energy E by the expressions k 2 — 2\i(E — 
t n )/h 2 and Vn = Z?Z% ii/{k n h 2 ). 

The partial wave functions ^^(r) satisfy the regularity 
condition at the origin 



*/i(0) = , 
while, as r — > oo, they have the asymptotics 



(10) 



M^(r) ~ SfifFtitH, h r) + u+( Vf ,k f r)A% . (11) 

Here the functions 

uf(ri n , k n r) = e Tia " [Gi(r) n , k n r) ± iFi{r\ n , k n r)\ 

are defined by using the regular (irregular) Coulomb 
wave function F e (Ge), and <r" = arg |r(l + I + irj n )} 
is the Coulomb phase shift in the n-th channel [23j. 
It should be noted that if the Coulomb interaction is 
not present in a channel then the Coulomb wave func- 
tions in the asymptotics (ITT1) should be replaced by the 
Riccati-Bessel ji and Riccati-Hankel hf functions [23[ 
since F e (0, k n r) = je(k n r) and uf(0, k n r) = hf(k n r). 

The quantities A% are the scattering amplitudes due to 
the short-range interaction V. Their dependence on I and 
E is assumed implicitly. The total scattering amplitudes 
are given by the sum 



A fi = A° r 5 



where 



A c 



f°fi 



A 



fi 



exp(2zcr/) — 1 



'It 



(12) 



is the partial Coulomb scattering amplitude. The partial 
wave cross sections are then determined by the amplitude 
through the standard expression 



7 Si 



4tt(2£+1)|%| 



(13) 



The total cross section corresponding to the reactive scat- 
tering transition i — > / is given by the sum over momenta 



tot 

a Si 



1=0 



(14) 



In the following discussion matrix notation is used. 
Here, the set of wave functions ^^-(r) are considered as 
a square matrix with the indices / and i running over all 
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values from 1 to M. All matrices are denoted by bold 
typeface. In this matrix notation, Eq. takes the form 



[H + L(r) + C(r)-fc 2 ]*(r) 

+V(r)*(r) = 



(15) 



Here k and Ho represent the diagonal matrices k n 5 nm 
and —d 2 /dr 2 5 nm , respectively. The centrifugal term 
L(r), the Coulomb interaction term C(r) and the short 
range potential coupling matrix V(r) are defined by the 
following matrix elements 1(1 + l)/r 2 5 nm , 2k n rj n /r S nm 
and V nm (r), respectively. The total Hamiltonian matrix 
is given by 

H = H + L(r) + C(r)+V(r). 

The regularity condition (ITU1) takes the matrix form 

*(0) = . (16) 

By introducing the diagonal matrix [u(r)] nm — 
u i (Vn,k n r)5 nm the boundary condition (|11|) can be 
rewritten as follows 



*(r) - F(r)+u{r)A s 



(17) 



where F(r) = [u(r)B - u*(r)]/(2i) and [D] nm = 
exp[2ier"]<5 nm . Due to the short-range interaction, V, 
the amplitude matrix A s is constructed from the partial- 
wave scattering amplitudes A s ^ with all possible values 
of indices / and i. The total scattering amplitude matrix 
is then given by the sum 



A = A 



A, 



(18) 



where the Coulomb scattering amplitude matrix (|12p is 
denoted as A c . 

Solving Eq. (|15H17[) by the exterior complex scaling 
technique requires a reformulation of the problem. In 
a similar approach to that employed in the one channel 
case 0, [H| the new incident wave which incorporates 
the long range tail of the interaction potential should 
be constructed on the first step. This approach will be 
described in the next two subsections. 



C. The solution to the problem of the long range 
diagonal part of the potential 

Here we consider the diagonal part of the equations 
fTS"]) with the interactions due to the long-range tail 



[Ho + L(r) + C R (r) - k 2 ] ^ R (r) = (19) 
with the regularity boundary condition 

ip R (0) = (20) 
and with the asymptotics as r — > oo 

i/> R (r) ~ F(r) + u(r)A R . (21) 



Solving the problem (|19II21[) is naturally reduced to the 
construction of the solutions to the individual equations 



1(1 + 1) 2k n rj n 



dr 2 



9(r-R)-k 2 n $?(r)=0, 
(22) 

where 0(t) is the Heaviside step function defined such 
that 0(t) = 0, t < and 9(t) = 1, t > 0. The diagonal 
Coulomb tail potential C R n (r) is represented in Eq. (|22|) 
by its explicit form 2k n i] n /r 9(i — R). The scattering so- 



lution to Eq. (|22l) can be constructed using the matching 
lur 

then ip R (r) takes the form 



ng the mai 

procedure that is described in detail in [J 01- 1$ r < R 

(23) 



^n( r ) =M(k n r)a R , 
or in matrix notation 

^ R (r)=l(r)a R 
and if r > R then ?p R (r) is given by 

ip R (r) = e ia ?Ft (n„, k n r) + uj(n n ,k n r)A 
or in matrix notation 

il> R (r) = F(r) +u(r)A R . 



R 

n ' 



(24) 



(25) 



(26) 



The matrices a R and A R in (|24l |26|) are diagonal. The 
values of diagonal elements a R and A R follow from the 
matching conditions at the point r — R and read 



(27) 



A 



exp [2iarg(a^)] - exp[2ia"] 
2i 



(28) 



where Wij(/, g) denotes the Wronskian f(r)g'(r) — 
f'(r)g(r) calculated at r = R for the functions / = 
u~l(r)n,k n r) and g = je(k n r). Finally, the diagonal ma- 
trix tp R (r) is defined by 



(29) 



W R (r)] nm = iP R (r)5 n 



Another kind of the solution to the equation ([19[) is 
defined by the asymptotic condition @ as r — > oo 



u R (r) ~ u(r). 



(30) 



This solution can be constructed by the same matching 
procedure just employed and this results in the following 
form of the components of u R (r) 



u n( r ) = uj(r] n ,k n r) 



(31) 



for r > R and 

u R (r) = h~(k n r)c R + h+(k n r)d R (32) 
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for r < R. The coefficients c R and d R are given by 

c R = W R (u+, hf)/W R (hi,hf) 
d R = W R (u+, hj)/W R (h+, hj), 



(33) 



a. R = a R 



where the Wronskians are computed for the functions 
u# (in, k n r) and hf(k n r) at the point r = R. The diag- 
onal matrix 



[u R (r)] nm = u R {r)5 n 



(34) 



provides the solution to the Eq. ([19)) . 

The solutions i^ R and u R allow us to construct the 
Green's function g R by the standard formula 



g R {ry) = k- 1 ii, H {r < )u H {r > ) 



(35) 



where r>(r<) = max(min)[r, r']. This is possible since 
the diagonal matrices ip R and u R commute. By con- 
struction this function obeys the equation 



[H + L(r) + C R (r) - k 2 ] g R (r,^ 



where the matrix has the form 

oo 

I- k- 1 J dr' u(r')V(r')^ R (r') . (42) 
R 

For r > R, the equation (|40)l reads 

r 

* fl (r) = xp R (r) - k^uir) J dr' tP R (r')V(r')H> R (r') 

R 

oo 

-fe-V fl (r) J dr'u(r')V{r')V R (r'). (43) 

r 

The asymptotic form of \I> fl (r) asr->oo can now easily 
be evaluated from the right hand side of this equation 
by neglecting the last term, since it goes to zero, and 
extending to infinity the upper limit of the integral in 
the second term. The result from this is given by 



= I5(r-r'), (36) ^> R {r) ^ ^ R {r) - k^u^) dr' tp R {r')V{r')<b R {r') 



where I denotes the unit matrix. 



D. The scattering problem for the entire potential 

tail 



Let us now consider the scattering problem for the en- 
tire tail potential matrix C R (r) + ~V R (r). The coupled 
Schrodinger equation in this case reads 



n 



(44) 



[H + L{r) + C R (r)-k 2 ]^ R (r) 
-V R {r)V R (r). 



(37) 



Regularity at r = and asymptotic as r — > oo, the 
boundary conditions take the form 



* fl (0) = , 

* fl (r) ~ F(r) + u{r)A R . 



(38) 
(39) 

The solution to Eqs. (|37M39|) are conveniently obtained 
from the solution of the Lippmann-Schwinger integral 
equation 

oo 

«r«( r ) = ij> R (r) - J dr'g R {r, r')V R (r')^> R {r'). (40) 



This equation uses both the solution xp R and the Green's 
function g R that are constructed in the preceding sub- 
section. Due to © these type of equations have a unique 
solution |24j . From this equation it follows that \& K (r), 
similar to xp R (r), takes a different functional form if 
r < R or r > R. 

Using (|23]l for r < R and the relevant representation 
of the Green's function g R (r,r') one arrives at the ex- 
pression 



From this formula the final asymptotic form (|39l) of 
*& R (r) as r — > oo can be obtained by using P5|). There- 
fore, the scattering amplitude A R is given by 

oo 

A R = A R - k- 1 J dr' [F(r') + u(r')A R ] V(r')H> R {r'). 

R 

(45) 

The Green's function G R (r, r') which satisfies the 
equation 

[H + L(r) + C R (r) - k 2 ] G R (r,r') 

+V R {r)G R (r, r') = IS(r - r') (46) 

can also be defined with the help of Lippmann-Schwinger 
equation 

oo 

G R (r, r') = g R (r, r') - J dr" g R (r, r")V(r")G R (r" , r'). 

R 

(47) 

This equation is also well defined due to ([8]) and has 
a unique solution. In the next subsection we use the 
asymptotics of the Green's function G R (r, r') in the spe- 
cial cases where r' < R and r ^> R. This asymptotics 
can again be evaluated from the right hand side of the 
equation (j47|) when r > R and r' < R 

G R (r,r') = k- 1 u(r)i/> R (r') 

r 

-k~ l u(r) J dr" il> R {r")Y{r")G R {r" ,r') 

R 



fe"V R (r) / dr"u(r")V(r")G R (r",r') (48) 



*«(r)=j(r)a«, 



(41) 



" „./~."wrf-*."\r~>Rf»." „'\ 
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by neglecting the last term, since it goes to zero as 
r — > oo, and by extending the upper limit of integration 
to infinity in the second term. This gives the following 
expression 



G fl (r,r') ~ fe- 1 u(r)* il (r / ), 



(49) 



where 



<& R (r') = xp R (r') 



dr" ip R (r")V(r")G R (r",r'). (50) 



By direct calculations one can verify that the transposed 
matrix v R (r) obeys the Lippmann-Schwinger equation 
(|40| for ^f R (r) and, therefore, due to the uniqueness of 
the solution of this equation the following equality holds 
true 



4f R (r) = * fl (r), 



(51) 



where T is the matrix transposition. The final form of 
the desired asymptotics where r' < R and r — > oo are 
obtained by taking into account the representation (|41[) 
and reads 



G R (r, r') ~ fe- 1 it(r)a iiT j(/). 



(52) 



E. The driven Schrodinger equation and exterior 
complex scaling. The integral and the local 
representations for the scattering amplitude. 

The solution *& R just obtained can now be considered 
as the incoming wave. By its construction, the action of 
the operator H — k 2 on this incoming wave has the form 



(H - k 2 )* R (r) = [Cfl(r) + Vfl(r)] * R (r). 



(53) 



If the "scattered" wave \& sc is introduced by the expres- 
sion 



*(r) = * R (r) + * sc (r), 



(54) 



then the Eq. ()15[) for ^(r) transforms into the inhomo- 
geneous equation for *& sc (r) 

(H - fc 2 )M> sc (r) = - [Cfl(r) + V H (r)] * fl (r). (55) 

This equation has two key properties which are very im- 
portant for application of the exterior complex scaling, 
i.e. the inhomogeneous term vanishes outside of the ra- 
dius R and the solution \I/ SC has purely outgoing asymp- 
totics. The final form for the driven equation formulation 
can be obtained from (|55|) by using the following observa- 
tion. The right hand side term in (|55|) has the following 
explicit form 

[C fl (r) + V fi (r)] * R (r) = [C R (r) + V R (r)}j(r)a R . 



By multiplying Eq. (|55[) by the inverse matrix (a. R ) 1 
from the right and then introducing the matrix 

*(r) = ^(rXa*)- 1 (56) 

the former equation transforms into 

(H-fe 2 )*(r) = -[C R (r) + V R (r)]j(r). (57) 



The boundary conditions follow from (fTo) . (fTT)) and 
(El 



*(0) =0 

<f>(r) ~ u(r)(A s - A R )(sl R )-\ 



r — >• oo. 



(58) 



Eqs. (|57l58p provide us with the final formulation for the 
driven Schrodinger equation. 

The integral representation for the scattering ampli- 
tude which results from the driven Schrodinger equation 
formulation is the last feature discussed in this subsec- 
tion. In order to derive this representation Eq. (|57|) 
should be recast into 

[H + L(r) + C R (r) + V R {r) - k 2 ] *(r) - 

- [C fl (r) + V fl (r)] [j(r) + *(r)] . (59) 

Using the Green's function G R (r, r') this equation can 
be rewritten in the integral form 

R 

*(r) = - J dr' G R (r,r')[C R (r') + V R (r')} 



jV) + *(0 . (60) 



The asymptotics of the solution <&(r) follow now from 
(|60| by taking into account the asymptotics (|52t of the 
Green's function G fl (r, r'). This gives 



*(r) - w(r)a fl J($) 



(61) 



J($) = -fc- 1 / dr' j(r') [C(r') + V(r')] [j(r') + *(r') 
o 

(62) 

By comparing (l58l) with ([6lj) and using the definition 
given in (|18p we obtain the final representation for the 
total scattering amplitude A 



A = A c +A H + a K ^ J($)a J 



(63) 



Thus if the matrices A R and a. R have been calculated 
then the driven equation formulation provides an alter- 
native to the original formulation (fT5)) , ([T^| and (fT7| . 
The important feature of this alternative formulation 
is that in order to calculate the scattering amplitude 
A one needs to know <&(r) only in the finite interval 
(0 < r < R). 
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The complex scaling application is based on the follow- 
ing arguments. The function <&(r) has purely outgoing 
wave asymptotics at infinity. When the complex scaling 
transformation is applied to the boundary problem given 
by ([57l [58]) , the scaled boundary problem will have the 
zero boundary conditions both at the origin and at infin- 
ity. The driving term in the driven Schrodinger equations 
(|57p does not diverge at large distances under the com- 
plex scaling transformation of the coordinate if, in the 
exterior complex scaling, a proper choice of the scaling 
point is made. Therefore, the necessary conditions for 
application of the exterior complex scaling to the bound- 
ary problem given by (|57l |58|) are fulfilled. The formal 
scheme of this application is as follows. If we denote the 
complex scaled function as $(r) = $>(g at Q(r)) then, for 
this function, we obtain the boundary problem 

' (H - fc 2 ) *(r) = - [C fl (r) + V K (r)] j(r) , 
* *(0) = , ( 64 ) 
k *(oo) = , 

where H represents the complex scaled Hamiltonian. The 
finite-range driving part remains unchanged after the 
complex scaling if R < Q. Furthermore, if R < Q then 
the complex scaling transformation does not change the 
value of the function 3>(r) in the region r < R and, as 
such, J(Q) = J($)- Thus, provided that we have solved 
the scaled problem (f6~4")l for <I>, the scattering amplitude 
matrix A can be computed from the representation 

A = A° + A R + h rT J($)a fl . (65) 

The matrices a R and A R defined by Eqs. flU]) and (143)) 
have to be determined in order to use (|65p . According 
to their definitions, these matrices can be calculated if 
we know the function ^f R (T) in the region r > R. How- 
ever, the numerical integration of the differential equa- 
tion (|37[) with an arbitrary potential V(r) and the point 
R is a problem of similar complexity to the initial scatter- 
ing problem (|15H17|) . Therefore, in practical calculations, 
we choose the point R to be large enough in order to 
assume that V fl (r) = 0. The validity of this assumption 
should be checked for each potential under investigation. 
The truncation of potentials decreasing faster than r~ 2 
at infinity does not lead to principal errors contrary to 
the truncation of the Coulomb potential. In section UlI Bl 
we will analyze how the truncation of the potential V 
affects the total cross section. 

If we set V R = 0, then Eqs. flUJ) and (05$ yield 

a R = a R , (66) 
A R = A R , (67) 

where the matrices a R and A R are diagonal. The inte- 
gral representation for the scattering amplitude (|65[) then 
transforms into 

A = A° + A R + a R J{$)a R . (68) 



Furthermore, the asymptotic relation (|6"Tj) becomes exact 
for r > R. Taking into account the fact that the function 
at the point r = R is not complex scaled, $>(R) = &{R), 
we conclude that 

A = A c + A R + u- 1 {R)${R)a R ; . (69) 

The last expression provides us with the local representa- 
tion for the scattering amplitude which is an alternative 
to the integral representation (|68j) . 

III. NUMERICAL APPROACH, RESULTS AND 
DISCUSSIONS 

The equation with zero boundary conditions (1641) to- 
gether with the two alternative representations for the 
scattering amplitude and can be directly im- 
plemented numerically. However, the derived equations 
can be slightly modified in order to obtain a more nu- 
merically stable implementation for large orbital momen- 
tum. Calculations with such momenta are necessary to 
achieve converged results for the total cross sections, see 
section UlIBl In the section Hi CI we described a method 
for constructing the solution to the diagonal part of the 
equation (|15p which incorporates the long range Coulomb 
interaction for r > R. Another approach is obtained 
through splitting the centrifugal term L(r) for each par- 
tial wave, 

L(r) = Lfl(r) + L H (r) (70) 

in the same way as is done for the Coulomb interaction. 
Eq. (fT9")) in this case transforms into 

[H + h R (r) + C R {r) ~ k 2 ] xP R (r) = 0. (71) 

functions hf(k n r) should be replaced by the trigonomet- 
ric functions sin(fc„r) and e ±lfe " r , respectively, in all of 
the formulas given in the section III CI Then the final 
equation (|64p transforms into 

J" (H - fc 2 ) *(r) = - [C«(r) + Mr) + V R (r)} sin(r) , 
\ *(0) = , ( 72 ) 
I *(«>) = , 

where [sin(r)]„ m = 5 nm sxa(k n r). The scattering ampli- 
tude matrix can then be calculated with the same local 
(|6"9")) and integral (f6"8"]l representations. However, the ma- 
trices and A R must be calculated using the expres- 
sions (|27|) and (|28|) with the above stated modifications. 

The boundary problem (|72| rather than (|64|) is used 
in all of the calculations reported in this section. For the 
sake of completeness, it should be noted that there exists 
a third possibility, where the unperturbed Hamiltonian 
also includes the Coulomb potential. Then the solution 
of the scattering problem with the long range diagonal 
part of the potential is explicitly given in terms of the 
Coulomb functions. This version has been explored e.g. 
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in papers pilOjj]. However, it is not clear how this tech- 
nique can be extended to systems consisting of three or 
more particles. Therefore, we do not focus on this option 
in this paper. 

We use the derived equations in order to study two sim- 
ple models. We consider first the short-range one-channel 
Noro- Taylor potential [25[ supplied with the Coulomb 
interaction. The second example is the more realistic 
two-channel model for the N 3+ + H -> N 2+ + H + re- 
action [l3|, This latter model is composed of the 
Coulomb interaction and molecular inverse power poten- 
tials. 

As the numerical method for the solution of Eq. (|72l) , 
we have chosen the FEM-DVR approach described 
in 28]. This approach can be considered as a finite el- 
ement method (FEM) with a special choice of the basis 
functions on each element, namely the polynomial Lo- 
batto shape functions (28|. With this choice, the matrix 
elements of local operators (i.e. potentials) are approxi- 
mately diagonal with respect to the basis functions. The 
error introduced by such approximation does not influ- 
ence the convergence rate of the FEM. The parameters 
of the numerical applications were chosen such that the 
numerical inaccuracies were negligible. 

Although the radius R and the exterior complex rota- 
tion radius Q are allowed to be different in the scheme 
described above, we have not found any advantages to 
keeping them distinct. Hence, in our calculations we 
choose Q = R. The preliminary calculations have also 
confirmed that the specific choice of the exterior com- 
plex scaling defined by Eqs. ([U and dl]) does not affect 
the results. Therefore, we have used the sharp exterior 
complex scaling [llj in the calculations. 

( a _ / r for r < R ,„ . 

9aM r ) ~ | R+(r-R) e «* for r>R' 

In contrast to applications of the complex scaling method 
to computing resonances, the choice of the scaling an- 
gle a is not limited here by any additional restrictions. 
Therefore, the angle has been chosen to be close to 90 
degrees to enable the fastest decay of the wave function 
at infinity. 

In order to employ the boundary conditions in Eq. (|72[) 
into the numerical scheme, we introduce the maximal ra- 
dius R m ax > R, where the second boundary condition of 
(|T2")) is implemented. The radius R ma x should be consid- 
erably larger than R, such that the wave function decays 
on the interval [R, R ma x]- As soon as R m ax is sufficiently 
large, no noticeably errors in the results are observed. 

A. Specific aspects of the calculations with the 
Coulomb potential 

Consider the scattering problem on the one-channel 
Noro- Taylor potential [1^] in the presence of the repulsive 
Coulomb interaction 

V{r) = I5r 2 e~ r , C(r) = 2/r. (74) 




^K) 20 30 40 50 

Radius R, a.u. 



FIG. 1. The partial wave cross sections <r°, a 5 , and a 10 (the 
dotted, dashed and solid lines, respectively), are plotted as a 
function of the radius R for the energy E = 3 a.u. 

Channel indices are not used since the model contains 
only one channel. The reduced mass was chosen to be 
(1 = 1. The total cross section for this potential with the 
Coulomb tail is infinite but one can analyze the partial 
wave cross sections. In the computations, we have used 
1000 finite elements with Lobatto polynomials of sixth 
degree. 

Let us first discuss the influence of the radius R on the 
partial wave cross sections. In Fig. [1] we plot the partial 
wave cross sections as a function of the radius R. As the 
short-range part of the potential (174]) decreases very fast, 
an accurate value for the cross section can be obtained 
already for rather small radii R starting from R — 20 a.u. 
These values depend on the energy E and on the short 
range part of the potential such that V s (R) <C E. The 
accuracy is not influenced by the value of the angular 
momentum £, at least for moderate values of I. 

The matrices a R and A R are defined through the ex- 
plicit expressions (|2"7T) and (|28p . The Coulomb wave func- 
tions are used in Eq. (|2T[) and in the local representation 
(|69[) . Keeping in mind the three-body generalization of 
our approach 10], we can check here whether we can 
completely avoid using the Coulomb functions in our cal- 
culations. 

The standard approach for this check is to use the 
asymptotic expansion of the Coulomb functions for large 
R. One can then show (see for example [23| for the one- 
channel case) that two first terms in the 1/R, R — > oo, 
expansion for u(R) and a R are given by 

[u{R)] mn ~ S mn (l + u n ) exp[i9 n ] , 
[a\ n ~ S mn (l - w„) exp [ir} n log(2k n R)} . (75) 

Here 9 n — k n R — ir£/2 — r\ n log(2fc n i?), and 

_ vn + tiiie + ^ + ti) 
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Energy E, a.u. 



FIG. 2. (Color online) The partial-wave cross section <r° as a 
function of the scattering energy E. The results for the exact 
boundary condition (|27[l (the solid line) , asymptotic boundary 
condition (|77[) (the dashed line), and the asymptotic bound- 
ary condition with the correction term (|78|) (the dotted line) 
are shown. The radius R — 100 a.u. 

ivl + Vn exp [2ik n R - ml] 
W " = 2%Ji • (?6) 

In the numerical scheme (fT2"j) , the similar asymptotics are 
modified to be equal to 

[u{R)]mn ~ S mn exp [i9 n ] , 
[a R ] mn ~ S mn exp [i (tm log(2fc„i?) + ni/2)] , (77) 

for the main terms of the asymptotics. The next order 
terms are given by 

[ aR ] mn ~ 6mn ^ ~ W «) CX P l0g(2/c„i?) + 7T^/2)] , 

ii£ + ?7„ exp [2ik n R] + i£(£ + 1) 
W " = 2fc^R ' (78) 

while it(i?) coincides with that in Eq. (|7rj|) . It is noted 
that the expressions (|77|> . ([75"|) . and ([75)1 are only valid 
when |u n | <C 1. Let us check how the derived asymp- 
totic representations influence the cross section calcula- 
tions. In Figs. [5] and [3] we plot the partial wave cross 
sections as a function of the scattering energy E. We 
compare the cross sections obtained with three different 
boundary conditions at the point r = R, namely the ex- 
act boundary condition (|2"7|) . the asymptotic boundary 
condition (|77[) . and the asymptotic boundary condition 
with the correction term (|75)) . For the chosen parameters 
and energy regions, the results for the integral and local 
amplitude representations are indistinguishable and, as 
such, we only plot the results for the integral represen- 
tation. For zero angular momentum, Fig. [5J all curves 
practically coincide. Some differences appear only for 
small energies, where the value for the asymptotic pa- 
rameter u given by Eq. (|76[) approaches |u| ~ 0.25 from 
below for E = 0.03 a.u. and gets even bigger for smaller 




Energy E, a.u. 



FIG. 3. (Color online) As in Fig. [2] but for the partial-wave 
cross section a 10 . 



energies. For the momentum I = 10, Fig.[3l the difference 
between the boundary condition (|2"T)) and the two other 
boundary conditions is more pronounced and clearly in- 
creases with decreasing energy. Conversely, the results 
for the exact boundary conditions (|27p and the asymp- 
totic boundary condition with the correction term (|78[) 
agree quite well, even starting from relatively small en- 
ergies E > 0.2 a.u. . On the contrary, the results for the 
simple asymptotics ([771) disagree with the correct results 
over the entire energy region shown. Thus, we conclude 
that the correction term introduced in (|78[) improves the 
cross section essentially when compared to the simple 
asymptotics (ITT)) . This also means that with the correc- 
tion term (|78[) we can, for a given energy, use a smaller 
value of R in order to reach the same accuracy. 

The accuracy of the correction depends on the param- 
eters u n defined in Eq. (|76[) . If the angular momentum 
is fixed, the accuracy is improved when kR increases. 
Conversely, for chosen scattering energy and radius R, 
the accuracy gets worse when the angular momentum 
increases. 

With respect to the differences between the inte- 
gral (|68)) and local (|69|) representations of the scattering 
amplitudes, they seem to depend on the chosen boundary 
conditions. For the exact boundary condition (f2~T[) both 
these representations give identical values, while they re- 
sult in different values for the asymptotic boundary con- 
ditions (|77l 1751) . This difference grows when kR gets 
smaller, as discussed in Ref. Q- 

Our approach is rigorous in treating the Coulomb in- 
teraction both in the one channel case 0] and in the 
multi channel case described in the present paper. The 
one channel scattering problem with a potential that de- 
creases at large separations faster than r~ 2 was treated 
in the paper [9(. The driven equation of the form (|64[) 
was used and then solved numerically. In the notations 
of this paper, the scattering amplitude of Ref. [9j is ex- 
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FIG. 4. The s-wave cross section a as a function of the 
scattering energy E. The cross-section is calculated using 
the exact integral representation (|68|) (the solid line) and the 
representation given in (|79l) (the dotted line). 



FIG. 5. The total cross section aX\ as a function of the energy 
E. The values for the radii R = 100 a.u. (the solid line), 
R = 70 a.u. (the dashed line), R = 50 a.u. (the dash-dotted 
line) and i? = 30 a.u. (the dotted line) are plotted. 



pressed as 

A = J(#) . (79) 

Let us compare this representation with the exact 
one (|68p. First of all, if the Coulomb interaction is not 
present in the system, then A c = and a R = I, A R = 0. 
Hence, the representations (|68|) and (|79|) are identical. If 
the Coulomb interaction is present, we can use the fact 
that the matrices a R and k are diagonal and so rewrite 
the exact matrix representation (|68l) in terms of the ma- 
trix elements as 

A fi = A c f 8 fl + AfS fi + afj(Q>) fi af . (80) 

The diagonal elements A^ and A R are equal to zero for 
the channels without the Coulomb interaction. There- 
fore, for all inelastic transitions and for elastic transitions 
within the non-Coulomb channels, we obtain 

A fi = af fi af . (81) 

If the Coulomb potential is absent in both the / and 
i channels, then a R — a R — 1, and there is no differ- 
ence between the two representations. If it is present 
in one or both of the inelastic channels then, according 
to the asymptotic representation (|75[) . \a R \ —> 1 when 
R — > oo. This means that the representation described 
in ([75]) gives the same cross section as the exact represen- 
tation (|rj8|) provided that R is sufficiently large. However, 
it should be noted that, as opposed to the cross section, 
the scattering amplitude is not calculated correctly even 
in the R co limit. 

In summary of this comparison, we can distinguish 
three different combinations of channels / and i for 
the Coulomb multi channel scattering problem. If the 



Coulomb interaction is absent in both channels (non- 
Coulomb elastic and inelastic channels), the results given 
by the representation shown in (|79|) are identical with the 
exact values defined with Eq. (|rj5|) . If the Coulomb inter- 
action is present in an inelastic channel, the partial cross 
sections calculated with Eq. (|79")l approach the correct 
values as R — > oo. For the Coulomb elastic channels, the 
representation (f79|) fails to give the correct answer. In 
order to illustrate the latter statement, we compare in 
Fig. H] the s-wave cross section calculated with both the 
exact integral representation (|68p and the scattering am- 
plitude (|79|) for the potential (|74|) . These cross sections 
completely disagree. 



B. The iV 3+ + H ->• iV 2+ + H+ reaction 

Here we study the two channel charge transfer 
N 3 +(ls 2 2s 2 ) + H(ls) ->■ NH 3+ ->■ N 2+ (ls 2 2s 2 3s) + H+. 
The matrix potential describing this reaction is parame- 
terized as 

, _ / 4000r- 8 - 20.25r- 4 0.5r 2 e - r \ 
V W =2/1 ( V 0.5r 2 e~'' -0.235 ) ' 

C(r) = 2 M (° . (82) 

The reduced mass was taken to be ijl = 1713.5 a.u. These 
parameters are taken from [27| while motivation for the 
choice of the model and its parameters are found in (2(| . 
In the numerical study of this system we have also used 
the FEM-DVR with the sixth degree Lobatto polynomi- 
als. The number of equidistant finite elements depends 
on the radius R. Their density has been chosen to be 12 
elements per atomic unit. 
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FIG. 6. (Color online)The total cross section cr*2* as the func- 



tion of the energy E. 
(the solid line) and R 



The values for the radii R = 50 a.u. 
= 30 a.u. (the dashed line) are plotted. 



In this discussion we will mainly concentrate on the 
development of the method to determine an appropriate 
choice of the only parameter in our approach, namely the 
cut-off radius R. In order to show the importance of this 
choice, we plot in Fig.[S]the total elastic cross section, <r*°* 
(IT4l) . for different values of the radius R. Analysis of these 
results show that although the form of the cross-sections 
are similar, calculations with relatively small values of R 
essentially underestimate the value of the cross section 
compared to those calculated with the largest value of 
R. Such a situation can easily result in large absolute 
errors in the computation of the the total cross section. 
The total rearrangement cross section, cr'^'i i s plotted in 
Fig. [6l In this channel, the total cross section converges 
already at the relatively small value of the cut off radius, 
R = 50 a.u. Hence we expect that the R value which 
guarantees the convergence depends on the behavior of 
the channel potentials at large distances. 

Let us therefore estimate the error which is introduced 
in the total cross section (fT4")l due to the choice of poten- 
tial cut-off. First of all, the cut-off results in a change 
of each partial cross section. In order to estimate this 
change, we should refer to the representation of the scat- 
tering amplitude (I63I) for the uncut potential. Assuming 
that the change of amplitude due to cutting the potential 
is much smaller than the amplitude itself, we can see that 
the main change comes from the last term in Eq. (I63I) . 
Therefore, for large R, we have 
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where 



oo 

W R = k- 1 J dr' u(r')V(r')V R (r'). (84) 



Using the asymptotic expansions for the functions u(r) 
and \& fi (r) for large r > R, for the matrix elements of 
the integral we obtain 



W 



R 



,„„ ~ . dr' k^e^ r '-^^V mn {r')s^{k n r' - n£/2). 
R 



(85) 

Let us assume that dV mn (r)/dr <§C V mn (r) for large r. 
For example, such an assumption is valid for the impor- 
tant type of potentials: inverse power potentials. Inte- 
grating ([55)1 by parts, we find for the main term of the 
W R n asymptotics at large R 

,i(k m ~k n )R. p i(km+k„)R—ini 



W R m — k^ 1 }/ 



for fc n k m , and 



w 



R 

mn 



iR 



J2ik m R-ml 



(86) 



(87) 



for k n = k m , and where the potential V mn (r) decreases 
as ~ C/r 13 at infinity. For the sake of simplicity, we 
omit here the indexes mn for the parameter ft. These 
expressions result in the components of the amplitude 
change 5 A 



[8A] fi 



(W R J nl (<P)a R a R + J fn (<S>)a R a R W R ) 




The corresponding change in the partial cross section (|13[) 
is given by 



(89) 



The representation (|88|) shows that, in general, the am- 
plitude error due to the potential cut-off depends on the 
amplitude in the various channels. For a specific sys- 
tem, however, some channel potentials can decrease much 
faster than others, resulting in a simpler description. For 
example, in the system the off-diagonal potentials 
decrease exponentially and this also means an exponen- 
tial decrease in W R and W R X . Thus, for the elastic non- 
Coulomb channel 1 —> 1, we find that a R — 1 and 

[6A] U w -2Jn($)W^. (90) 

Taking into account the asymptotics (1871) . we obtain for 
the relative error of the partial cross section 
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for large R. In this case, the relative error decreases as 
r~> R~ 3 . The absolute error of the total cross section is 
estimated as 
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FIG. 7. The partial wave cross section o~n as the function of 
the momentum I for the energy E — 2 • 10~ a.u. is plotted 
for different values of R: R = 500 a.u. (the solid line), R — 
100 a.u. (the dashed line), R — 75 a.u. (the dash-dotted line), 
R — 50 a.u. (the dotted line) and R = 30 a.u. (the dashed 
line). 



The absolute error also decreases with the same rate ~ 
i?~ 3 for large R and its value can be easily controlled in 
the calculations. 

Comparing the results plotted in Fig. [5] and given by 
the estimation (j92|) . one can see that Eq. (|92|) essentially 
underestimates the error. In order to clarify this differ- 
ence, let us discuss the dependence of the partial- wave 
cross section a\ x on the orbital momentum I. The corre- 
sponding values are plotted in Fig. [7] for a few values of 
the cut-off radii. The calculation for R = 500 a.u. can 
be considered, in the range of the figure, as a reflection of 
the exact results. Analysis of the results shows that there 
are three distinct regions of ^-dependence. For small val- 
ues of £, this dependence is not regular. For medium 
values of £ one can see a region of linear dependence that 
results in the inverse power behavior of the cross section 
as a function I. Finally, for large values of £, the par- 
tial wave cross sections decrease very fast. So while the 
results for the values of £ in the medium region can be 
considered as good approximations to the correct cross 
sections, the calculation for the large £ fails completely. 

In order to analyze Fig. [7J we should recall the am- 
plitude definition (|62)) . For a given cut-off radius R and 
large £, such that k n R < £, one can use the asymptotics of 
the Riccati-Bessel function with respect to its index [23[ : 
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(93) 



This function decays extremely fast, so does the integral 
«/(<!>). However, the exact scattering amplitude for the 
uncut potential has different behavior. In order to find 
its asymptotics, we notice that the expression (|83[) does 
not describe the leading asymptotic term as the condition 



8 A <C A is no longer satisfied because of the fast decrease 
of J(&). Conversely, the leading asymptotic term is now 
the integral term in the representation (l4"5j) for A R : 



A 1 



dr' k- 1 3t(k m r')V mn (T')Uk n r'). (94) 



For the 1 — > 1 channel, the potential Vu consists of a 
superposition of the inverse power potentials 1/r 13 {82}. 
The integral with standing instead of V\\ can be 

evaluated explicitly [2{|. This leads to the following 
asymptotics of the partial wave cross section (flU} at 
the large values of the momentum £ in the case of the 
inverse power potential 
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(95) 



The constant Cp is known explicitly and depends only on 
the potential. As the Vu component of the potential (|B"2"|) 
decreases as ~ 1/r 4 , this results in the crfj ~ l/£ 5 asymp- 
totic behavior of the partial wave cross section. This 
behavior explains the linear intermediate region in the 
partial wave cross sections in Fig. [JJ especially so for the 
semi-exact results with R = 500 a.u. 

The numerical results, however, are computed with 
the cut-off potential so the amplitude (|94| is not taken 
into account. So when £ becomes sufficiently large that 
the behavior of the function jg (r) approaches the asymp- 
totics the behavior of the calculated partial cross 
section changes. They vanish very fast when compared 
to the correct values, and this can be seen in Fig.[7J When 
summing up the partial wave cross sections to obtain the 
total cross sections, this might give rise to the erroneous 
impression that the total cross section is already con- 
verged. This implies that such incorrect behavior can 
influence the accuracy of the calculated total cross sec- 
tion. Furthermore, we cannot improve the accuracy by 
increasing the number of partial wave cross sections taken 
into account since the higher partial-wave cross sections 
are not correctly calculated. 

Let us estimate the sum of discarded partial wave cross 
sections for the inverse power potential 1 /r' 3 . We denote 
by £r ~ ekiR/2 the value of £ when the asymptotics (j9"3")l 
should be taken into account. The correction to the total 
cross section Acr'°' can then be estimated by using (|95|) 
as 



At 
<7j 



tot 
1 



E 



'ii 



Cpkt 



2/3-6 



oo 

t=t R 



,3-2/3 



2(0-2) 



(96) 

For a given inaccuracy in the total cross section Aer^f*, 
the latter estimation can be considered as the equation 
for the minimal value of R which guarantees the re- 
quested accuracy. According to the representation (|63p . 
the amplitude inaccuracies ([55)1 and are to be 

summed in order to obtain the total estimation. This 
implies that the double sum of partial wave cross sec- 
tions and (|95p gives the estimate for the inaccuracy 
in the cross section. 
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We would like to stress that the estimates (|92l |96| are 
specific not only to our approach. Any numerical method 
which cuts the potential in one way or another, experi- 
ences these errors. This implies that estimates similar 
to Eqs. (|92l [96]) should be adopted in order to guarantee 
the accuracy of the numerical calculations. Finally, we 
have compared our two-body partial wave and total cross 
section results with those described in ref. [l3| using a 
log-derivative method. We find an agreement within the 
accuracy discussed. We also find that the needed com- 
putational resources are comparable. 

IV. SUMMARY 

Inspired by the work of Nuttall and Cohen [8] and 
Rescigno et al 0] we have developed a theory which en- 
ables the calculation of two-body multi channel charged 
particle scattering. This theory is formulated with the 
aim of being able to generalize from two-body to three- 
body charged particle multi channel quantum scattering, 
as briefly outlined in [Icj . 

The entire potential is split into a core and a tail po- 
tential. The scattering problem in the first step is solved 
for the diagonal, analytically solvable part of the tail po- 
tential. The solutions of this problem are then used to 
construct the corresponding Green's function. These are 
then used to derive the Lippmann-Schwinger equation, 
the solution of which is the wave function correspond- 
ing to the entire tail potential, including the off-diagonal 
elements. This resulting wave function is then used as 
an incident wave when formulating a driven Schrodinger 
equation, which has the desired scattering function as 
a solution. Using exterior complex scaling, this scatter- 
ing wave function can be obtained from a boundary value 
problem (|64|) with zero boundary conditions at origin and 



at infinity. The scattering amplitude is then defined by 
Eq. (USD. 

The theoretical results are supported by the numerical 
realization using the FEM-DVR technique. The theory 
is illustrated by application to both a one-channel as well 
as a two-channel problem which both include Coulomb 
interaction. Our formulation for the problem with the 
Coulomb interaction is theoretically as well as numeri- 
cally compared to the one expired by Rescigno et al [9|. 
For multi channel scattering our analysis shows that the 
approach for the amplitude extraction of Ref Q gives 
the correct results for the cross sections in non- Coulomb 
channels and asymptotically correct results for inelastic 
Coulomb channels. However, the representation for both 
the cross section and the scattering amplitudes which 
follows from the Rescigno et al formulation cannot be 
directly generalized for the Coulomb elastic channels. 

For the practical implementation of our approach, we 
have introduced a cut-off radius for the short range po- 
tentials. We have therefore estimated errors in the total 
cross section due to this cut off. For the important class 
of inverse power potentials, we have found simple esti- 
mations for the minimal cut-off radius which guarantees 
a desired accuracy. 
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